Approach of a class of discontinuous dynamical systems of 
fractional order: existence of the solutions 
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In this letter we are concerned with the possibility to approach the existence of solutions to 
a class of discontinuous dynamical systems of fractional order. In this purpose, the underlying 
initial value problem is transformed into a fractional set-valued problem. Next, the Cellina's 
Theorem is applied leading to a single-valued continuous initial value problem of fractional 
order. The existence of solutions is assured by a Peano like theorem for ordinary differential 
equations of fractional order. 
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Let us consider a general class of autonomous discontinuous dynamical systems of fractional order 
modeled by the following Initial Value Problem (IVP) 



Dlx = f{x) := g(x) + A s(x) 



(0) 



(fc = 0,l,...,M -1), t€/=[0,oo) 



(1) 



where g : W 1 —> W 1 is a function continuous with respect to the state variable, A = ((H^) a real constant 
matrix and s is a piecewise continuous function given by 



/ sgn( 



s(x) 



\sgn(x n ) 



q G M+ and D% is considered in this letter as being the most utilized differential operator, the Caputo 
operator of order q with starting point 0, i.e. (see e.g. |Podlubny[[l999l ) 



Dlu{i) 



( t _ r )M-<?-i n (M)( r ) dr . 



r : (0,oo) 



1 « 

r(M -q) Jo 

is the known Euler's Gamma function and [•] denotes the ceiling function that rounds 



up to the next integer. Thus, Di 9 ^ is the conventional differential operator of order \q] G N. 



Remark 1. According to the standard mathematical theory Stefan et al., 1993, §42], we are forced to give 
the initial conditions for the IVP ([!]) using fractional derivatives of the function /, or these values are 
frequently not available. Also, it may not even be clear what their physical meaning is. Therefore, using 
the Caputo's suggested way, the initial conditions may be specified in the classical way, as in IVP 0. 
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For the sake of simplicity, we restrict ourselves to the case important for the applications: q £ (0, 1) 
(however the considerations in this paper can be generalized to arbitrary positive q). Therefore, we deal 
with the following form of IVP ([I]) 

D%x = f (x) := g (x) + A s{x) , x(0) = x , t G I = [0, oo), (2) 

where we have to specify just one condition since it is easily seen that the number of initial conditions that 
one needs to specify in order to obtain a unique solution is \q] =1. 

Remark 2. The IVP (fiT) are enough general to include the great generality of systems: for q = 1 the systems 
modeled by the IVP 111) are the known Filippov systems Filippov 1988 , while for non integer values of q 
and A = O nxn , the IVP ([I]) models dynamical systems of fractional order. 

The main result of this letter is the following theorem 

Theorem 1. The |I]) admits at least one solution. 

For a better readability of the letter, the proof of this theorem shall be given in several steps. 



1. Filippov regularization of the right-hand side 

The existence and uniqueness of solutions to discontinuous IVPs are essential for discontinuous dy- 
namical systems, because due to the right-hand discontinuity, classical solutions of IVP might not even 
exist. To provide the existence, it is necessary to modify the right-hand side of IVP ([2]). For discontinuous 
vector fields, existence and uniqueness of solutions is not guaranteed in general, no matter what notion of 
solution is chosen. Also, the classical notion of solution for ordinary differential equations is too restrictive 
when considering discontinuous vector fields. A possibility to compass this difficulty is to extend the notion 
of differential equation to differential inclusion, problem solved by Filippov using a generalized concept of 
solution. Thus, the single valued discontinuous IVP is shifted to the following set-valued one 

D\x £ F(x), x(0) = xq, for almost all t S I, (3) 

where F : K n => M n is a set-valued vector function defined on the set of all subsets of M. n . One of simplest 
definition of F is the following convex form (implicitly used in most introductory references) 

F( x ) = f] p| convf((x + eB)\M), (4) 

e>0 fi (M)=0 

where, M is the set of discontinuity points of /, B the unit ball in W 1 , \x the Lebesgue measure and 

conv the closed convex hull. At the points where the function / is continuous, F{x) will consist of one 
point, which is the value of / at this point, i.e. F(x) = {/(x)}. At the discontinuity points, the set F(x) 
is given by Q. Therefore, F(x) is the convex hull of values of / (x*), x* 6 M, ignoring the behavior on 
null sets. For example, the Filippov regularization applied to the unidimensional sign function leads to the 
set-valued function 

r {-1} x<o, 

Sgn{x) = I [-1, 1] x = 0, 
({+1} x>0. 

and the right-hand side of the IVP ^ becomes 

F(x) := g(x) + A S(x), with S(x) = (Sgn(x{), Sgn(x n )) T . (5) 

In order to justify the use of the Filippov regularization in physical systems, e in Q must be small enough, 
so that the motion of the physical system can be arbitrarily close to a certain solution of the differential 
inclusion. 



2. Continuous approximation of the right-hand side 

Next, X and Y denote metric spaces (e.g. W 1 as in almost real applications). 
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Definition 1. A selection of a given set- valued function F : X Y is a function h : X — > Y satisfying 

Vx G X, h(x) G F(x). 

Definition 2. A set-valued function F : X — > Y is called upper semicontinuous (u.s.c.) at x G X if for 
any neighborhood V of F(x), there exists a neighborhood U of x such that .F(x) C V for all x G V. F is 
u.s.c. on X if it is u.s.c. on every point of X . 

For practical reasons, it is convenient to characterize a set-valued map F : X ==>- Y by its graph 

Graph(F) = {(x,y) e X x Y \ y £ F(x)}. 



Proposition 1. The set-valued function F defined by is u.s.c. with nonempty closed and convex values. 



Proof. The proof can be found e.g. in Aubin & Cellina, 1984, P- 102]. 



Remark 3. Due to the symmetric interpretation of a set-valued map as a graph (see e.g. [Aubin 
Frankowskal [1990] ) we shall say that a set-valued map satisfies a property if and only if its graph sat- 



isfies it. For instance, a set-valued map is said to be convex if and only if its graph is a convex set. 

The following known theorem (Cellina's Theorem or "Approximative Selection Theorem") will be a main 
tool used in the proof of Theorem [l| 



Theorem 2. 'Aubin & Cellina, \1984\\Aubin & Frankowska , 19901 Let F : X — > Y be upper u.s.c. set- 



valued with Y a Banach space. If the values of F are nonempty and convex then, for every e > 0, there 
exists a locally Lipschity function f e : X — > Y such that 

Graph(f e ) C Graph{F) + eB. 

The proof of Theorem [2] is constructive (see the mentioned reference) in that it provides a method to 
explicitly construct the selection f £ with e parameter. This is an important advantage which allows to find 
practical approximations in real applications. 

3. Existence of the solutions for fractional equations 

Let us consider the following IVP of fractional order with q G (0, 1) 

Dlx = f{x), x(0) =x , (6) 

with x G M n . The following existence result corresponds to the classical Peano existence theorem for first 
order equations 



Theorem 3. Kai, 2010^ Assume f in IVP (|) is continuous and bounded. Then there exists a solution to 
IVP®. 



Proof, see Kai 2010 Corollary 6.4 p. 92]. 



Next we can give the proof for Theorem [T] 

Proof. Applying Filippov regularization, the right hand side of ([!]) transforms into the set- valued function 
F given by ([5]). F is a convex u.s.c. (Proposition [T]) and non-empty valued function (Remark [3]). Therefore, 
Theorem [2] can be used and the IVP ([T]) becomes a continuous IVP of fractional order to which Theorem 
[3] applies and the proof is complete. ■ 

Remark 4- i) Using the Fillipov regularization, the single valued IVP may be considered as transforming 
into a set-valued IVP of fractional order. However, for practical purposes, fractional single- valued IVPs are 
more accessible by numerical point of view, since there are several ways to approximate their solutions. 



ii) The uniqueness for general fractional equations is treated in Kai 2010, §6.2] the underlying theorem 
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corresponding to the well-known Picard- Lindelof Theorem for equations of integer order and is based on 
Lipschitz continuity. However, for our case of systems, the uniqueness is checked since from Theorem [2] 
under some a priori boundedness, the approximate selection is locally Lipschitz continuous. 



(7) 



For example, let us consider the following fractional discontinuous variant of the Chua's system 

Dlx\ = —2.57xi + 9 x 2 + 3.87 sgn (xi), 
Dlx 2 = x\ - x 2 + x 3 , 
Dix 3 = -px 2 , 

where p is a real parameter. The graph of the first component of the right-hand side of ([7]), the scalar 
function f\ (x\, x 2 ) = — 2.57xi + 9 x 2 + 3.87 sgn (xi), is plotted in Fig. [jja). The convex hull of /i (0, x 2 ) is 
the dashed region and represents the graph of the set- valued function i<i(0, x 2 ). The underlying set- valued 
IVP is 

Dlxi G -2.57xi + 9 x 2 + 3.87 Sgn (xi), 

Dlx 2 = xi - x 2 + x 3 , (8) 
Dlx 3 = -px 2 . 

A possible approximation (selection) of F\ in the neighborhood of (0, x 2 ), for x 2 £ R, can be for ex ample 
a smooth cubic polynomial function h (xi,x 2 ) = ax\ + bx\ + cx 1 + 9x 2 |Danca &: Codreanu 



2002 



where 

a, b, c and d have to be determined from the continuity and differentiability conditions in xi = ±e (Fig. 
0b). 

Next, the system can be either numerically integrated using e.g. the fractional Adams-Bashforth-Moulton 



e.g 



et al. 2002 . 



method discussed in |Diethelm et ah 2002 , or approximated using some frequency-domain method ( see 
Charef et al, 1992| ). More approximations ways for fractional operators can be found e.g. in |Vinagre 



Graph{FA0,.K 7 )) 




(a) 




Fig. 1. a) Graph of F\ (0,2:2); b) Approximation of F\. 



REFERENCES 5 



References 

Aubin, J. -P. & Cellina, A. [1984] Differential inclusions set-valued maps and viability theory (Springer, 
Berlin). 

Aubin, J. -P. & Frankowska, H. [1990] Set-valued analysis (Birkhauser, Boston). 

A. Charef, A. Sun, H. H., Tsao, Y. Y. & Onaral, B. [1992] "Fractal system as represented by singularity 
function," IEEE Trans. Automat. Contr. 37(9) 1465-1470. 

Danca, M.-F. & Codreanu, S. [2002] "On a possible approximation of discontinuous dynamical systems," 
Chaos, Solitons & Fractals 13 (4), 681-691. 

Diethelm, K., Ford, N. J., & Freed, A.D. [2002] "A predictor-corrector approach for the numerical solution 
of fractional differential equations," Nonlinear Dynamics 29(1-4), 3-22. 

Filippov, A. F. [1988] Differential equations with discontinuous righthand sides (Kluwer Academic Pub- 
lishers, Dordrecht). 

Diethlem K. [2010] The Analysis of Fractional Differential Equations; An Application- Oriented Exposition 

Using Differential Operators of Caputo Type (Springer- Verlag, Berlin). 
Podlubny, I. [1999] Fractional Differential Equations (Academic Press, San Diego ) 

Samko, S. G., Kilbas A. A. &: Marichev, O. I., [1993] Fractional integrals and derivatives: Theory and 

applications (Gordon and Breach, Yverdon). 
Vinagrc, B. M., Podlubny, I. Hernandez, A. & Feliu, V. [2000] "Some approximations of fractional order 

operators used in control theory and applications," Fractional Calculus & Applied Analysis 3(3), 

231-248. 



